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Abstract 

We propose a physically motivated parametrization for the unpolarized generalized parton dis- 
tributions, H and E, valid at both zero and non-zero values of the skewness variable, £. Our 
approach follows a previous detailed study of the £ = case where H and E were determined 
using constraints from simultaneous fits of the experimental data on both the nucleon elastic form 
factors and the deep inelastic structure functions in the non singlet sector. Additional constraints 
at C 7^ are provided by lattice calculations of the higher moments of generalized parton distribu- 
tions. We illustrate a method for extracting generalized parton distributions from lattice moments 
based on a reconstruction using sets of orthogonal polynomials. The inclusion in our fit of data 
on Deeply Virtual Compton Scattering is also discussed. Our method provides a step towards an 
extraction of generalized distributions based on a global fit of the available data within the given 
set of constraints. 
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I. INTRODUCTION 



A completely new perspective was given to the study of hadronic structure, with the 
observation that deep inelastic exclusive experiments, such as Deeply Virtual Compton 
Scattering (DVCS), and hard Exclusive Meson Production (EMP), could in principle allow 
one to access spatial configurations of partons besides the relatively well known longitudinal 
momentum fraction ones 1|. Because of the specific role played by the partons' spatial coor- 
dinates, one can also envisage progress towards the solution of the proton spin crisis through 
studies of the orbital angular momentum contribution of different partonic components |2|. 
Information from exclusive deep inelastic scattering is coded in terms of Generalized Parton 
Distributions (GPDs) representing the soft matrix elements in the process. Since they were 
first introduced 0, 3, 4], much progress has been made in determining their general prop- 
erties in terms of the relevant kinematical variables, including the longitudinal momentum 
fraction taken by the struck quark, X, the four-momentum transfer defining the scale of the 
deep inelastic process, Q 2 = —q 2 , the four- momentum transfer squared between the initial 
and final proton states, t = A 2 , and the longitudinal momentum transfer fraction of the 
initial proton momentum, the "skewness", £, or £, (for reviews see 

At present, a central issue is the definition of a quantitative, reliable approach beyond 
the construction of GPDs from specific models and/or particular limiting cases, that can 
incorporate incoming experimental data in a variety of ranges of Q 2 and A = (£, t), as they 
gradually become available. One hopes eventually to obtain a similar sophistication level 
as for the Parton Distribution Functions (PDFs) extracted from inclusive Deep Inelastic 
Scattering (DIS). We address this issue in the present paper. 

The matching between measured quantities and perturbative QCD (PQCD) based pre- 
dictions for DVCS/EMP proceeds, in principle, similarly to the inclusive case, with a few 
important caveats due to the fact that GPDs describe amplitudes and are therefore much 
more elusive observables from the practical point of view of experimental measurements, 
than the PDFs. In DIS the cross sections at large enough Q 2 and s (the total invariant 
mass), where a partonic picture is expected to hold, measure directly the Bjorken x depen- 
dent structure functions; x is interpreted as the momentum fraction carried by the partons, 
modulo target mass corrections. A widely used approach is "global fitting" whereby PDFs 
are obtained from experimental data by assuming a parametric functional form at a given 
scale, which is: i) evolved to the scale of the data; ii) convoluted with the appropriate per- 
turbatively calculated coefficient functions; Hi) compared to the data, thus determining the 
parameters. 

In DVCS/EMP similar factorization theorems as for the inclusive case exist. However, 
the phase space that has to be covered in order to extract GPDs has two extra dimensions, ( 
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and t, and the cross sections are written in terms of convolutions of GPDs over the variable 
X. Experiments delivering sufficiently accurate data have just begun. The comparison 
with experiment and the formulation of parametrizations necessarily encompasses, therefore, 
other strategies using additional constraints, other than from a direct comparison with the 
data. 

Here we explore the constraints on the extraction of GPDs provided by a combination of 
experimental data on nucleon form factors, PDFs, and lattice calculations of Mellin moments 
with n > 2. The latter, parametrized in terms of Generalized Form Factors (GFFs), were 
calculated by both the QCDSF 

BEL BEE Q and LHPC Q, 0, Q collaborations for 



both the unpolarized and polarized cases up to n = 3, 1 therefore allowing to access the 
skewness dependence of GPDs. We implemented results from Refs. [lo|, [Ll| using two flavors 
of (9(a)-improved dynamical fermions. Lattice QCD calculations are increasingly improving: 
new results can, and will be added to our determination as they become available. 

In order to estimate GPDs by knowing a finite set of Mellin moments we adopt a strat- 
egy similar to the ones originally developed in QCD based analyses of DIS data where a 
mathematically rigorous method for inverting moments was devised [3]. The early QCD 
studies were motivated by the simplifications that Mellin moments offered in treating PQCD 
evolution. In our case, moments provided by lattice QCD represent, in addition an impor- 
tant source of information on GPDs. We envisage a future scenario in which the calculated 
Mellin moments will be embedded in a "global" analysis including determinations of GPDs 
from experiment, and a consistent treatment of PQCD evolution. 

The problem of the construction of a function from its moments is mathematically mean- 
ingful so long as one establishes the necessary and sufficient conditions on the moments that 
a unique solution exists. A point-by-point reconstruction is therefore unattainable, however 
practical methods can be applied to a number of physics problems including QCD based 



analyses of the DIS structure function data [15| (see also [16j for a review of early analyses), 



so long as one aims at obtaining F 2 in a finite interval around XBj- It should be noticed, 
in fact, that experimental measurements always provide discrete "binnings" of kinematical 
variables. As shown in Refs. 15| the Bernstein polynomials are ideal for reproducing the deep 
inelastic structure functions in that they are zero at the endpoints, they are normalized to 
one, and they are peaked in different regions within the interval xbj £ [0,1]. Because of 
the latter property the Bernstein polynomials allow one to emphasize the behavior of the 
structure function at given specific regions of xsj, while suppressing the others. It was found 
that n > 8 moments were necessary to give a fully quantitative description of the behavior 



1 In principle, lattice calculations can extend to n = 4. This has not been done so far because of statistical 
noise. 
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of F2(xBj,Q 2 )- At present, only n = 3 Mellin moments are available from lattice QCD cal- 
culations, therefore one cannot reach a similar level of accuracy as in the inclusive case. In 
this paper we show that we can, however, reduce the kinematical region where information 
comes strictly from this type of reconstruction u sing; Mellin moments solely to X < (. The 



parametrization we obtained in the £ = case 17| can in fact be extrapolated smoothly 
to X > £ (the so-called DGLAP 18| region) by a simple kinematical shift in the quarks 
ongitudinal momenta, being based on the same partonic picture. The X < ( - ERBL 
19| - region does not lend itself however to a clearcut interpretation in terms of quark and 
gluon structure. It is in fact becoming increasingly manifest that in order to extract and 
interpret the amplitudes from a variety of hard exclusive experiments including DVCS in 
the ERBL region, it is fundamental to understand and connect the partonic structure of 

|. The 
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21 



t-channel exchanges to previous information from Regge phenomenology 
contribution of a Regge-motivated term was already considered in quantitative studies at 
0, 



C 



23. 



24] . following a similar logics as in the early ansatze for PDF parametrizations 



(see e.g. [25|] ) . Regge type contributions are expected to play a much more important role 
for GPDs at ( 7^ 0, specifically in the X < ( domain, due to the dominance of the scattering 
from a qq pair emerging from the initial nucleon in these processes (FigQ~b). 

The large theoretical uncertainty in the X < ( domain motivated us to perform our anal- 
ysis concentrating on a determination of GPDs in this region, dictated only by experimental 
data and lattice moments, and therefore in principle model independent. The price that one 
pays for maintaining model independence is, however, in this case, and with the present ac- 
curacy of lattice QCD determinations, a loss in predictivity in the point-by-point functional 
dependence of the GPDs in the ERBL region. Our results however provide points from 
ab initio calculations, with given error bars in a region that cannot be extracted directly 
from experiment. These can be used in turn for guidance to constrain the shape of possible 
functional forms of GPDs, a question addressed in a forthcoming manuscript 36j. Finally, 
improvements in the accuracy of the extracted points will be made as more accurate lattice 
QCD results become available. 

An important consideration should also be given to the analytical properties of GPDs at 
X = Q. This point marks the transition between the partonic and mesonic interpretations 
of GPDs, respectively, and it acquires a particular significance only in the presence of soft 
matrix elements describing the non-perturbative structure of hadrons (in e.g. a QED based 
description it no longer signifies a specific breakpoint). Simple consistency arguments related 
to factorization ^] require that continuity is observed at the breakpoint. Recent precise 
measurements from Jefferson Lab indeed find non-vanishing values of GPDs of the measured 
Single Spin Asymmetry (SSA) at X = ( 26j. The continuity condition constitutes an 



additional constraint in our approach in addition to the ones provided by the experimental 



4 



data and the lattice results. 

It should be remarked that our results are specifically for the flavor Non-Singlet (NS) 
sector. Perturbative evolution is therefore slower than for the singlet and gluon cases, and 



27 



28]. 



it can be implemented straightforwardly similarly to the PDF case 

The main thrust of the present paper is to present an extension of the analysis of Ref.[l7| 
using available information at £ ^ 0. Our method can be seen as a first step in a "bottom up" 
approach as opposed to other models presently available {e.g. the dual model of Ref. [2l|, and 
the Double Distribution based model of Ref.[29|). Our aim is to make use of all information 
available both experimentally, and from ab initio calculations, and to study how to construct 
the functional dependence of the GPDs consequentially. As in similar phenomenological 



studies (see e.g. [23|, |24|, [30|) a simple and adequately flexible functional form given by the 
spectator diquark model is used in the partonic-process-dominated DGLAP region, with 
the goal of eventually easing into a global analysis. As for the ERBL region, our scope is 
to define what constraints are available, that come directly from ab-initio calculations. A 
similar exploratory study has not been conducted so far to our knowledge. The study of a 
functional form in X in the ERBL region that is consistent with the constraint found in our 
first quantitative study is the subject of a forthcoming manuscript 36( where the real part 
of the DVCS interference term will be presented. As a first step, we preferred limiting the 
present manuscript to providing a conistent method that uses all information obtained from 
experimental constraints plus lattice QCD results. Future extensions to the sea quarks and 
gluon distributions will also help disentangle the intricate interplay between perturbative 
evolution and shape of the initial distributions. 

Our paper is organized as follows: In Section II we define the ingredients of our physically 
motivated parametrization in the ( > case, and introduce lattice moments with a prescrip- 
tion to account for their chiral extrapolation. In Section III we describe the extraction from 
lattice moments using Bernstein polynomials; Our results are compared to the recent data 



from Jefferson Lab 



In Section IV we draw conclusions and give an outlook. 



II. CONSTRAINTS ON GENERALIZED PARTON DISTRIBUTIONS FROM 
LATTICE MOMENTS 



A. Generalities 



GPDs can be extracted most cleanly from Deeply Virtual Compton Scattering (DVCS) 
at four-momentum transfer, Q 2 , in the few/multi-GeV region (Fig{T]). For unpolarized scat- 
tering one has two independent functions: H, and E, from the vector (7^) and tensor (07^) 
interactions, respectively. These are parametrized in terms of three kinematical invariants, 
besides the initial photon's virtuality, Q 2 . In this paper we adopt the following set: (£, X, t), 
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where ( = Q 2 /2(Pq) is the longitudinal momentum transfer between the initial and final 
protons (( ~ xbj in the asymptotic limit, with Bjorken XBj = Q 2 /2Mi>, v being the energy 
transfer associated to 7*), X = {kq)/(Pq) is the momentum fraction relative to the initial 
proton carried by the struck parton with momentum k, t = A 2 , is the four-momentum 
transfer squared. In principle, the observables for DVCS and similar exclusive processes 
are functions of only two independent invariants: t and s (s = Q 2 (l/xBj — 1) + M 2 ). 2 
A partonic picture is however more conveniently described in terms of Light Cone (LC) 
momentum projections such as ( and X, despite the fact that X is not directly observable, 
and it appears in the amplitude as an integration variable HQ. The need to deal with 
a more complicated phase space, in addition to the fact that DVCS interferes coherently 
with the Bethe-Heitler (BH) process, are in essence the reasons why it is more challenging 
to extract GPDs from experiment, wherefore guidance from phenomenologically motivated 
parametrizations becomes important. 



Initial experiments were conducted at HERMES 3l|. However, more recently, the high 



precision measurements from Jefferson Lab 26j have provided both the real and imaginary 



parts of the Bethe-Heitler BH-DVCS interference term at Q 2 ~ 2 GeV 2 , and XBj = 0.36, 
showing the dominance of the twist-2 contribution. This, in turn, can be written as j^ : 

C x = Fi(t)W(C, t) + ^G M (t)H((,t) - J^F 2 (t)£((,t). (1) 

In Eq.(TjQ), T = {H, £, H}, q = u,d, s, define the generalized Compton Form Factors (CFFs) 
for the unpolarized and polarized case, respectively; ^1(2) are the Dirac (Pauli) form factors, 
Gm = F\ + F2 (for details on the photon electro-production cross section and the harmonic 
decomposition of its (BH) 2 , (DVCS) 2 , and BH-DVCS interference terms we refer to Refs.jsl, 



3*. 

Several definitions have been given connec ting the CFFs to GPDs (see the review in [5|]). 



We follow the notation and definitions of Ref. 33J] according to which the leading order CFFs 
can be written as: 



HC, t) = -it £ e 2 q [F»(C, C, t) - F«(-C, C, *)] + 



Notice that Eq. (T5]) is completely analogous to the definition in the symmetric scheme of j^J , 
except for the choice of kinematics. For practical calculations it is usfeul to consider the 



2 Throughout the paper we use the "asymmetric" notation with skewness parameter £ = Q 2 /2(Pq). This 
is related at leading order to the skewness parameter. £. in "symmetric" notation by: £ = C/(2 — C)- I n 
this representation X = {x + £)/(! + £)■ 
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separation into valence and sea quarks GPD components. Specifically, GPDs are related 
to scattering amplitudes, at variance with the PDFs, and their quark/ ant i- quark content 
is a function of the longitudinal momentum transferor skewness, (. At non-zero skewness 



different definitions were given in the literature (e.g. |28l. |34J|). The essence of the argument 
is that for GPDs three separate kinematical regions define the scattering process: 

i) £ < X < 1, where the initial quark in Fig{T]has momentum fraction X and the final one 

x-C\ 

ii) 1 — ( < X < where one scatters an initial anti-quark with momentum fraction —(X — 
> 0, into a final one with —X > 0. This corresponds to a ^-channel exchange (the lower 
limit of integration in Eq.([2]) is defined accordingly); 

Hi) < X < C, where a quark- antiquark pair emerges from the nucleon and undergoes the 
electromagnetic interaction. Barring contributions from very low X (X < 0.01) dominated 
by anti-quarks in the initial state, in this region the quark carries momentum X > 0, and 
the anti-quark — (X — Q > 0. 

In summary, the amplitudes for scattering from either a quark or an anti-quark are 
defined in different kinematical regions. As for the identification of the valence and sea 
quarks contributions, analogously to the PDFs 35], (by omitting the t-dependence that is 
not relevant here), one has for the forward (£ = 0) case 



F q (X,0) = q(X) 
F q (-X,0) = -q(X), 



(3a) 
(3b) 



q{X), and q{X) being the quark and anti-quark parton distributions from DIS. Therefore, 
the valence forward GPD, F y (X,0), is 



F^(X, 0) = F q (X, 0) + F g (-X, 0) = q(X) - q(X), 



(4) 



and, following Ref. 35J, the sea quarks distribution coincides with the anti-quarks one, 



F q s (X,0)=q(X). AtC>0, 



F q (X,() = F q v (X,C)+F g s (X,0 X>0 
F q {X,C) = if(X,C) A>0, 



(5a) 
(5b) 



with — F q (— X, () = F q (X, () for — 1 + ( < X < 0. Once the physical meaning of the various 
regions in X is clarified, it is therefore a matter of choice whether or not to represent the 
different valence and sea quarks components of GPDs on the positive side of the X axis. 
Eqs.f l3ll4l[5l) are in line with Ref. 34|, where definitions where given in the context of DVCS 
from nuclei, that avoid the appearance of spurious symmetries in the X dependence of the off 
forward distributions [28| . Notice the consequences of the valence and sea quarks separation 
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on the structure of Eq. ([2|) : while the first term includes a sum over all q and q distributions, 
the second term is given by a principal value integral where the valence and sea quarks 
appear in the different domains, i), ii), Hi), explained above. In particular, since the valence 
contribution is zero at X < 0, this implies that the 1/X term in the intergral corresponds 
to a singularity that cannot be handled by the integration (a similar condition appears for 
x = — £ in the symmetric scheme). As a consequence, either the valence contribution to 
the CFF cannot be calculated alone, or the condition above dictates the dependence of the 
valence contribution to the GPD at X — > 0, which should be steep enough to counter the 
1/X type singularity. We explore this point in a forthcoming manuscript 36j. Finally, as 
observed in the forward case, F//(X, () ^ Fj?(X, (), in general. This is an aspect to be still 
explored e.g. for scattering from strange quarks in the exclusive channels. 

We also notice that a contamination from the pure DVCS term might be present in 
the extraction of the coefficient that might be large, although weighted by a kinematical 
factor of about 1% at Jlab kinematics. We give our estimate of this term in future work. 
Models such as do not seem to accurately reproduce the data on either the imaginary 
or real part of Eq.([T]). Although future experiments are planned that will allow one to 
directly determine whether such discrepancies could be also due to the presence of a pure 
DVCS contribution so far disregarded in the extraction of C, it is at present important 
to provide GPD parametrizations using constraints from both inclusive deep inelastic and 
elastic scattering data, and ab initio lattice calculations. It is our aim to explore, in what 
follows, the usage of such constraints towards a parametrization built from a bottom-up 
approach. 



B. Method Description 



In a previous publication 17|] we presented our parametrization of unpolarized GPDs 
in the flavor Non Singlet (NS) sector based on a diquark spectator model improved by a 
Regge term at low X. Parameters were given for skweness ( = 0, by using constraints from 
both form factors and PDFs. The spectator model is ideal in this context because despite 
its simplicity, it has proven to be sufficiently flexible to describe (and predict) the main 
features of a number of distribution and fragmentation functions in the intermediate and 
large X regions, as well as the unintegrated PDFs 
into a quark carrying a LC momentum fraction X - 



37 



391 ] . At X > C, the proton splits 
transverse momentum k i , and 



a spectator system with 1 — X, = k^-/P + , and — k 



k+/P 

. After undergoing the electromagnetic 



interaction, the final quark with momentum fraction X—(, and the spectator system coalesce 
into an outgoing proton (all particles are moving forward). The partonic configurations are 
therefore equivalent to the ( = case, modulo a kinematical shift. The parametrization's 
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form is therefore (see 171]): 



H{X,C,t) = G{X,(,t)R{X,(,t), 

where G(X, £, t) is the diquark model component, and R(XX,t) the Regge-based term. 
The parametrization differs from the one in 17J only by a shift to non-zero skewness, in the 
kinematical variables (a similar form is obtained for E(X,£,t)). 



Notice the difference with other parametrizations reproducing the £ = behavior |23|, [24] , 
where the forward contribution is given directly by existing PDF parametrizations, q(X). 
The latter are factored out from so-called "profile" of the parton distributions, namely: 



H{X,0,t) = q{X)exp[-tf{X)]. 

While in this case the forward limit, H(X, 0, 0) = q{X), is automatically met by introducing 
explicitly q(X), in our case, this needs to be enforced non trivially. 



In other words, in Ref. 



17] we had to simultaneously obtain new parametrizations for 



both the GPDs and their forward limit, the PDFs at a low initial scale. This effort granted 
us, however, the possibility to automatically extend our parametrization to £ ^ (X > Q 
by a simple variation in the kinematics. More details will be given in Section IIHI 

We subsequently analyzed the additional, ^-dependent constraints given by the higher 



re NS combinations: H u d = H u — H d , and 
1 corresponding to the nucleon 
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11], n 



moments of GPDs. The n — 1, 2, 3 moments of t 
E u-d = E u_ E d are avai i a ble from lattice QCD 

form factors. We use such constraints within a reconstruction procedure for the GPDs from 
their moments, using Bernstein polynomials. Such a procedure allows us to extract values 
for H u ~ d and E u ~ d at three values of X, defined as Xk,2, k = 0, 1, 2, for each given ( and t. 
The Bernstein polynomials are used as weighting functions for H u ~ d and E u ~ d , emphasizing 



the regions around Xk,2 Ha, lla]- As we will show in detail in Section [TTT1 moments of GPDs 
calculated using Bernstein polynomials provide the values of H u ~ d and E u ~ d , defined as 
H^. d and E^. d , at the points X^p, k — 0, 1,2. Numerical results including the errors and 
dispersions in X, respectively, are also provided in Section [TTTl 

Finally, we observe that the lattice results on GPDs moments have to be chirally extrap- 
olated. For this purpose we extended to the n = 2, 3 moments a simple ansatze proposed in 



41] for the nucleon form factors. The procedure is illustrated in Section III CI 

Because only the first three moments are known to date, by performing our analysis in 
the whole region of X e [0, 1] one obtains large theoretical uncertainties in the evaluated X 
and H, E values. 

We therefore determine the behavior of the GPDs at X > (, by keeping the same param- 
eter values obtained using the form factors and PDFs constraints in Ref. 171 ] . We checked 
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as explained later on that the curves at X > ( are consistent with the values obtained with 
the reconstruction procedure in the whole X interval. 

The Bernstein polynomials reconstruction is subsequently used in a reduced region, the 
"unknown" X < (, ERBL 19] region, where one expects a definite departure from the 
partonic type description of the DGLAP region. As we shall see in Section 1111} reducing the 
size of the interval in X provides an advantage in the GPDs extraction. 

Below we summarize all the information that was used in our analysis: 

• Accurate experimental information on the nucleon form factors. 



Parametrizations of the NS part of PDFs (e.g. the Alekhin set from Ref. 40|) repro- 
ducing the DIS data. 

Lattice results for the n = 1, 2, 3 moments at — t < 2 GeV 2 with the following features: 
— Only non-singlet contributions, H u ~ d , E u ~ d are provided (see Section Til C\} . 



The contribution from the term Co 



see 



zero, consistently with lattice calculations 



q s.fll2lll3p below), is set equal to 



12 



13|). 



— The dependence on ( is extracted from the term oc A 32 (see Eq. ( 118ft below). 

It should be noticed that experimental data on DVCS can in principle be systematically im- 
plemented in our extraction and will be considered in a forthcoming manuscript. Predictions 
for the data in Ref. 



261 ] are given in Section IIHI 



We conclude by underlining once more that our goal is to provide a practical method for 
extracting GPDs deriving from techniques that have been well tested in DIS. The proposed 
method uses information from both experimental data and lattice results, and it is exact in 
that for each kinematics, or ((,t) values, it provides the values of H(X,(,t) and E(X,(,t) 
in at given X values, with calculable theoretical error. With the proposed approach we also 
wish to provide an alternative based on a bottom-up type of analysis, at variance with top- 
down models proposed so far including t 
the Mellin-Barnes integral representation 



re Double Distribution (DD) hypothesis 29], on 



421 ] . and /or on the dual representation of Ref. 21]. 



Our approach, similarly to what found for DIS 15j] provides a different perspective on the 
problem of a formal extraction of GPDs from their moments that tends to be more subject 
to uncontrolled numerical ambiguities such as the one in e.g. the oscillating term inherent 
in the integral defining the continuation to complex n. 



C. Extrapolation of Mellin Moments from Lattice QCD 

The Mellin moments of GPDs are most clearly described in the symmetric frame of 
where the relevant kinematical variables are x = (k' + + k + )/(P' + + P + ) = ((A — £/2)/(l + 
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C/2) and f = -2A + /(P'+ + P+) = (/(2 - (). The x moments of GPDs are defined as 

H*(£,t) = j\xx n - l H«{x,t„t) (6) 

K&t) = J dxx^E'faU), (7) 

where q = u,d,s and we disregarded strange quarks contributions. H% and E% represent the 
form factors of local twist two operators. From the Lorentz structure of these form factors, 
one obtains the following polynomiality relations in £, made explicit by using the Gordon 
decomposition: 



n-l 
2 



H*{£, t) = + mod (^ 2)e n C«(t) (8) 

i=0 

n-l 

Kit,*) = E^(*)^-mod(n,2)r^(t). (9) 



i=0 



A^ 2i (t), B^ 2i (t), and C^ 2i (t) are the Generalized Form Factors (GFFs). The latter were 
calculated on the lattice by both the LHPC and the QCDSF groups for the unpolarized, 
polarized and transversely polarized cases 0, S, @, 0, Q- Since the GFFs correspond 
to off-diagonal matrix elements of twist-two operators, it was necessary to introduce a new 
technique implementing al 
the same continuum GFFs 



if (4) cubic group operators and index combinations producing 



121 ] in order to provide statistically accurate lattice measurements 



from an overdetermined set of equations. Therefore, lattice calculations extend only up to 

n < 3. We list below the GFFs expressions for n — 1, 2, 3 using the notation of 0, \\2 \. 

For n = 1 one obtains the nucleon Dirac and Pauli form factors: 

i i 

Hf = A\ Q (t)=F!(t) = J dxH*(x,£,t) = J -^H^XX.t) (10) 

-i -l+C 2 

1 1 

E\ = B\,(t) = Flit) = J dxE*(x,Z 1 t)= J ^E^Ct). (11) 

-i -l+C 2 

For n = 2 one has: 

Hi = AUt)+(-^) 2 Cl(t) (12) 

El = Bl (t) - ("^) 2 C q 2 (t) (13) 
In the limit t — > 0, Eqs. (110|llj) give the baryon number, B q , and anomalous magnetic 
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moment, K q : 



S?o(0) 



B n 



dXq(X) 



(14) 
(15) 



Eq. Q12p gives the LC momentum fraction carried by quark q: 



i 

A1 (0) = (x), = y rfXXg(X) (16) 
o 

Furthermore, Eqs. (ll2|ll3l) are related through the angular momentum sum rule 



Al (0) + BU0) = 2J«, 



(17) 



where J q = L q + S* 9 is the total - orbital plus intrinsic - angular momentum carried by the 
quark q. 

Finally, the n = 3 moments are given by: 



m 

El 



AUt) + eAUt) 



esiM 



(18) 
(19) 



In our calculation we use lattice calculations for the unpolarized GFFs for n = 1, 2, 3 ob- 
tained by the QCDSF collaboration using two flavors of (9(a)-improved dynamical fermions 
for several values of t in the interval < t < 5 GeV 2 , and covering a range of pion mass 
values, run > 500 Me V 2 0, [ill- Similarly to previous evaluations 12] the GFFs for both 
H and E display a dipole type behavior for all three n values, the value of the dipole mass 
increasing with n. 

A straightforward implementation of present lattice calculations in realistic parametriza- 
tions is hampered by the rather large discrepancy with the experimental data, associated 
to the large values of the pion mass, m n used in Refs. 1 121 . Il3| . Early extrapolations 
used a linear approximation ll| that although improving the comparison with experiment, 
is not sufficient to grant the accuracy that is necessary for a quantitative parametrization. 
Extrapolations using Chiral Perturbation Theory (ChPT) of the GPD moments up to n = 2 
are currently being addressed using different methods: Heavy Baryon ChPT (HBChPT), in 



Re 



46 



3 



44] . self-consistently improved ChPT [45|], and finite-range regulators techniques 
While most results have bean focused on the values at t — 0, which are important 
for the determination of the nucleon orbital angular momentum, equations for t ^ were 
given explicitely in 44J. Whether this can be applied to the large masses where most 



12 



lattice results are provided, is still a matter of intense debate that is beyond the scope of 
this paper. In addition, the number of parameters of the calculation is rather large to enable 
a really precise evaluation of the moments at t 7^ 0. Nevertheless evaluations providing the 
t" term with a ~ 25% error have been possible within the most recent set of lattice results 



481 ] . On the other side, the finite-regulator technique results by allowing for an extension 
of the range in pion mass considered, will conisderably reduce the uncertainty on the n = 2 
moment determination, While a number of new extrapolation methods are currently be- 
ing explored we adopted a simple ansatz for the extrapolation of the dipole masses for the 



Dirac and Pauli form factors [4l| which: i) uses the connection between the dipole mass and 
the nucleons radius; ii) introduces a modification of the non analytic terms i the standard 
chiral extrapolation that suppresses the contribution of chiral loops at large m n . Despite 



its simplicity the conjecture of [4l[ well reproduces the trend of lattice results at large 
while satisfying the main physical criteria i) and ii). The resulting values for the "physi- 
cal" nucleon dipole masses are both sensibly different than using a linear extrapolation, and 
much closer to the experimental data. As a consequence, the proton and neutron magnetic 
form factors, as well as the proton electric form factors can be reproduced quite accurately, 
whereas non-negligible discrepancies were found only in the case of the more elusive neutron 
form factor. Based on the good agreement with experiment obtained using this extrapo- 



lation method, we extended it to the recent lattice data on n = 2, 3 GPD moments [10]. 
We reiterate that the choice of this simple method in spite of recent ChPT developments 
is justified in the context of this work by the fact that on one side we need all available 
higher moments, up to n = 3, and these are not accessible with available methods; secondly, 
the recent flurry of activities has been focused on the t = region, while we are primarily 
interested in t ^ where the uncertainties in the various ChPT-based determinations be- 
come more important. Finally, the goal of our work is to suggest an approach to determine 
GPDs using all available constraints from both experiment and theory. While hoping that 
both present and future efforts will provide us with increasingly precise results, and with 
a method to determine larger n moments, the role of this work is limited to implementing 
them within our suggested analysis as a theoretical input with a precisely determined un- 
certainty, independently from the merit of any specific approach used to obtain them in the 
first place. A detailed study of different sources of theoretical uncertainties is beyond the 
scope of this work, and will be performed in a forthcoming paper. 

The chiral behavior of the form factors is exhibited by the isovector and isoscalar contri- 
butions: 

^M(E) = ^ V M(E) ~ ^M(E) (20) 
Gm(E) = ^ P M(E) + ^M(E)' (21) 
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where G P ^ E ^ are the magnetic (electric) form factors for the proton (neutron). 

Lattice results are given only for the isovector components, related to to the corresponding 
form factors as 

G V M = (Ff - F?) + {Ft - F?) (22) 
dX(H u -H d )+ [ dX(E u -E d ) 



and 

G V E = (Ff - F?) - t{F* - F 2 ") (23) 



J dX{H u -H d )-r J dX(E u -E d ), 



with r = —t/AM 2 > 0. From Eqs.fl22) and ([23]) one obtains 

iv i nv 



dX(H u -H d ) = E (24) 

E^ d = [ dX(E u -E d ) = G u~ G z . (25) 

J l + T 



( H M(E))n(t) = ~ 7~2 ' ( 26 ) 



For £ — it was found that the n = 2, 3 moments from Ref. [si, [l2] display a dipole 
behavior: 

Ji M(£)' 

where in this case: 

Hl~ d = j dXX n {H u -H d ) = t ( H m^ + ( H e)" (27) 

E u n - d = j dXX n {E u -E d ) = i^Mkz^k ^ (28) 

where the l.h.s. quantities are obtained from the lattice moments calculations, whereas 
(■^M(E))n an< ^ i^M(E))n are amenable to chiral extrapolations. 

We surmise that for n — 2, 3, similarly to the n = 1 case, a relation of inverse propor- 
tionality exists between the dipole mass and a given hadronic length scale, similar to the 
following relation between the magnetic and electric dipole masses, K\j and A E , and the 
corresponding mean squared radii exists, namely: 

<^V) = 7T^, (29) 

and that a chiral behavior for this hadronic scale can in principle be formulated. Lacking 
an exact expression for the coefficients in such an extrapolation, we performed a fit using 
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n 


(A}0 2 (GeV 2 ) 


(A^) 2 (GeV 2 ) 


1 


0.457 ± 0.048 


0.576 ± 0.060 


2 


0.704 ±0.163 


1.371 ±0.238 


3 


1.80 ±0.57 


1.25 ±0.43 



TABLE I: Values of the dipole masses squared, (A]0 2 and (A]^) 2 , for the vector electric and vector 
magnetic moments, for n = 1,2,3, as obtained from the fit of the lattice results described in the 
text. 



the following form: 



(A 



M{E)> 



i2(i±«r (i V: 



/3f (E) ±7nln 



(30) 



In Eq.( [30i) otn^ E \ fin and 7„ are parameters to be fitted to the lattice results; fi = 0.4—0.5 
GeV is a "cut-off" parameter for the pion mass, above which chiral loops are suppressed. 

Eq. (l30|) differs from the form factor expression from Ref. 41[ in that we let the coefficient 
of the non linear term, j n , be a parameter of t 
expression for this term was obtained in Ref. 
being the axial coupling constant, and f n 
Eqs.(j2 



re equation whereas in the n = 1 case, an exact 
3, namely: 7 i = g 2 A M N /(8<Kf^ v ), g A = 1.26 



93 MeV being the pion decay constant. 3 
(128]) were first fitted using^a dipole form for A^ d (t) = A u nG {t) - A d (t) and 

attice results. For n = 1 



B^ d (t) = B% (t) - B* (t), n = 1, 2, 3 |gj|, for all provided sets of 
the lattice results are consistent with what was obtained in Ref. 



( A £(M)) 2 ' are shown 



4l| . The values obtained at 
the physical pion mass are given in Table [H The dipole masses at n = 2 
in Fi g |2] along with our results for the fits using Eq. (l30~l) . with n = 0.4 GeV, in accordance 
with [4l| . 

FigIS] shows the moments A^ d (t), for n = 1,2,3, obtained with the values from Table 
[H For n = 1 the extrapolated lattice results are compared to the parametrization from 
Ref. 17|] that accurately reproduces the experimental data for the electric and magnetic 



form factors. Although there is a discrepancy with the parametrization/experimental data, 
its size is consistent with the previous evaluation from Ref. 41|. The error band on the 
lattice evaluations is a result of the fit error on the dipole mass. Finally, we also extracted 
the value for A 32 (t) = 0.0302 ± 0.0103 (Eq.(HHD), by performing a linear fit of the lattice 
results. 

We also note that similar results can be obtained for both E and the polarized GPDs, 
and will be presented in a forthcoming paper [3^] . In the next Section, in order to illustrate 



3 Note that a term oc arctan(/i/7r), present in 4l| was dropped for simplicity from the isovector magnetic 
term. The effect of disregarding this term for the first moment was evaluated to be within the errors from 
the lattice calculation. 



15 



our method, we concentrate on the function H, which displays the highest accuracy, being 
constrained in the forward limit. 



III. RECONSTRUCTION FROM BERNSTEIN POLYNOMIALS 

With a finite number of moments in hand, one can apply to GPDs a reconstruction 



method following the approach initially used for PDFs in Ref.[15|. The procedure, that 
uses only moments with integer n, introduces a complete set of positive definite polynomi- 
als, the Bernstein polynomials, from which one constructs weighted averages of the GPDs, 
around definite values of X. As pointed out in [lsl ]. because of inherent uncertainties in 
the experimental binnings, one does not need a point-wise description of the structure func- 
tions/GPDs, but averages of these quantities over ranges of X, thus making the method 
both mathematically sound and appealing from the practical point of view. The method is 
in fact an alternative to a reconstruction of the GPDs from their inverse Mellin transforms 
that however requires in principle the knowledge of the behavior with n of the Mellin mo- 
ments. Anyhow, such approaches were used for a quantitative reconstruction of PDFs from 



the deep inelastic scattering data 50j; it was more recently used in NNLO analyses in 51], 



and they were also suggested in order to evaluate GPDs at ( = [ll|. In a more general 
case of ^-dependent GPDs, merging information from both experimental data and lattice, 
and taking into account PQCD evolution, clearly involves a more complex phase space than 
for PDF studies, and makes the implementation of the inverse Mellin transform method a 
prohibitive task. 

We consider therefore: 

i 

Hk,n(C t) = J H(X, C, t) b Kn {X)dX k = 0, ...n, (31) 
o 

where the weight functions are the Bernstein polynomials: 

~y-k /i Y"\n—k 

h, n (X) = i [ } , (32) 

/ x k (i - xy- k dx 

o 

that, by definition, are functions that are concentrated on restricted ranges within the X G 
[0, 1] interval, specifically around the values: 

i 

X k ,n = J h, n {X) XdX = £±i (33) 



1(3 



with a width given by the dispersion Afc >n [15|, ll6j 4 - 

(k + l)(n-k + l) 



^k,n 



V2 _ V 2 
k,n ^k,n 



(n + 2) 2 (n + 3) 

By writing bk, n in Eq. fl3TT) with a binomial expansion, one obtains: 

H f , n _ (^ + l)' y (-1)' rj 



(34) 



k\ 



(35) 



1=0 



where Hi + k are the Mellin moments: 

Hi+k+i 



H(X,(,t) X l+k dX 



(36) 



We tested the range of validity of the proposed reconstruction approach by using an already 
known function, namely our parametrization for H(X,(,t) at C = [l7]. In Figf4]we show 
H(X,0,t) along with a reconstruction using the first eight (n = 7) and three (n — 2) 
Bernstein moments, respectively. From the figure it is clear that it is sufficient to consider 
eight moments to accomplish an accurate description of H u ~ d , at X > 0.1, and for all 
values of t. With only three moments, slightly larger discrepancies arise. It should be also 
noticed that the Reconstruction using Bernstein polynomials is always less accurate at small 
X, a feature already noticed in [3]. Such a discrepancy is however small with respect to 
current error bars on the lattice moments determinations. Moreover, if necessary, i.e. as 
more accurate lattice results and data become available, it can be improved upon either by 
im plem enting; higher order corrections 15| . or by using modified Bernstein functions as e.g. 



m 



52|. 



In Figures [5] and [6] we use the three available lattice moments (Section lH CI) to reconstruct 
H u ~ d . In this case, for each value of ( and t, one obtains the following values of X: 

X 2 = 0.25, X 12 = 0.5, X 22 = 0.75, 
and weighted average values of H: 



H 2(X 02 ) = 3A w -6A 20 + 3 
H 12 (X 12 A 
H 22 {X 22 ) = 3 A 30 



A 



30 



2C 



2-C 



30+ ( 


' 2( 


v 2-C 


-AH 


to 




A 32 





.4 



32 



(37a) 
(37b) 
(37c) 



Notice that the dispersion in Eq. (|M]) is different from the definition used in Ref. 52 1 
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In Figj5] we show the results for the weighted averages, with their given width calculated 
according to Eq. fl34"l) . The full curve in the DGLAP region was obtained by extending to 



C 7^ the parametrization from Set I of Ref. [17]] as follows: 

H T (X, C, t) = (X, C, t) R{{X, C, t) (38a) 
E 7 (X,C,t) = KG X M Ix (X,(,t)Rl(X,(,t) (38b) 



where G x Ml is given by a spectator model including both scalar and axial-vector components 
for the final diquark: 

Here k and k' are the initial and final quark momenta respectively (FigfjTJi), D(X, C,kj_) = 
k 2 — m 2 , D(X, (, k' ± ) = k' 2 — m 2 , k' ± = kj_ — (1 — X)/(l — ()A, m being the struck quark 
mass, k the quark's anomalous magnetic moment, A = P — P' being the four-momentum 
transfer, and: 

k 2 = XM 2 - ^M 2 X - (40) 



2 



k' 2 = ^4m 2 - ^M 2 X - k ± - i^A ) ±4, (41) 



with M, the proton mass, and Mx the (flavor-dependent) diquark mass (we suppress the 
flavor indices for simplicity). Furthermore, <p(k 2 ,\) defines the vertex functions in both 
the scalar and axial- vector cases 12J . The normalization factor includes the nucleon-quark- 
diquark coupling, and it is set to TV = 1 GeV 6 . (X,(,t) reduces to the form given in 
Ref. 17| in the £ — > case. 

The functions R{^ 2 ) were introduced in 17| in order to simultaneously account for: i) the 



low X behavior through an extra Regge motivated, term, X~ a , which is fundamental for 
obtaining the correct baryon number and momentum sum rules; ii) the small t behavior, 
while preserving the asymptotic behavior of the form factors. For ( = they read: 

R[ {2) =X- aI -^- x ? i{2)t (42) 



Notice that the resulting parameters for R\r 2 ) found in Ref. 17J are sensibly different for H 
and E, since at variance with the forward case governed by G T M where no specific constraint 
exists for E, these two functions are impacted on differently by the Dirac and Pauli form 
factors, respectively. Moreover, the fit parameters are not numerically similar to the ones 
from Regge phenomenology since the diquark term defined in Eqs.f l38a]) . fl38b[) keeps on 
being important at low Xand t thus interfering with the Regge part of the parametrization. 
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The overall behavior of H is however consistent with the parameter values from Regge 



phenomenology (see discussion in 17] and 23J). Finally, it should be remarked that Eqs.fl38 



present an approximate treatment of the spin structure, and that improved spectator models 



such as the ones recently worked out in Ref. 38j] should be adopted for a full consideration 



of all GPDs, and of their spatial d.o.f. interpretation. This is true in our case also for the 
function E where we chose a rather crude approximation, justified in our case by the goal of 
our paper that largely makes use of moments of GPDs. It should be kept in mind that this 
approximation is expected to eventually break down at very large X, as it already can be 
seen in Ref. 17j in reproducing sensible values for the transverse displacements s q . A full 
treatment of the spin structure of GPDs is however beyond the scope of this paper and it 
will be considered elsewhere. 

For C 7^ 0, R[(2) is a l so modified with respect to the ( = case through a ^-dependent 
kinematical shift in t resulting from the fact that 



t = + (1 - e 2 )D 2 ^ -t mm + (1 - e 2 )D 2 (43) 



where the Fourier conjugate to the impact parameter is now D = P'/(l — £) — P/(l + £) 



and f = C/(2 - C) (D reduces to A at ( = |53j), yielding: 



R{(XX,t) = x-^l^-xf^it+t^) (44) 

More details on this transformation are given in Appendix |A] From Figj5] one can notice a 
very good agreement between the Bernstein reconstructed values, and the parametrization 
in the DGLAP region obtained by perturbatively evolving at Leading Order (LO) Eqs. fl38al) . 
( I38bl) to the scale Q 2 = 4 GeV 2 , where the lattice results where evaluated. We also show 
the missing area (shaded area in the figure) obtained through the polynomiality condition, 
by ensuring that the sum of the areas in the ERBL and DGLAP regions, respectively, gives 
the form factor value at the given value of t. 

The Bernstein polynomials based reconstruction was subsequently repeated in the ERBL 
region only. This was accomplished by assuming the validity of the parametrization in the 
DGLAP region, and by calculating reduced moments in the ERBL region through: 

(H n ) x< = Hi = H n - (H n ) XH , (45) 

where (H n ) x>( . is given by: 



(Hn) 



1 

J H I (X,(,t)X n ~ 1 dX n = 1,2,3 (46) 
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with H ! (X, C, t) given by Eq. fl38al) . while H n are the lattice moments. The procedure defined 
in Eqs.( l37a]) was repeated for the reduced moments, yielding: 



#02(^02) = ^<|34 C 2 -64 C + 3 



^12(^12) = ^|e4 c-6 



H22{(X 2 2) = —z < 3A 30 + 



^30 



2( 



4 C 

^30 



2( 



2C 



2-C 



.4 



32 



2-C 



.4 



32 



2-C 



.4 



32 f ? 



(47a) 
(47b) 
(47c) 



Eqs.f l47p were obtained by redefining the Bernstein polynomials in the region X £ [0, Q as: 

X k ^ _ x yn-k 



bk,n(XX)= , 

/ x k (c - xy- k dx 

o 

the dispersion (Eq. (l34|) . being evaluated as: 



(48) 



(49) 



Results using the reduced moments are presented in FigJHl along with our evaluation in the 
DGLAP region using Eq. (j38ap . The missing areas, which value is given by = H^, are 
also highlighted in the figure. The theoretical errors displayed in the figure are from two 
sources: the lattice errors, and the fit errors (see Figj3]), added in quadrature. One can see 
that the effect of reducing the interval over which the weighted areas are taken, produces 
much larger errors on the moments, due to the presence of powers of 1/C in Eqs. (I47D . 
The determination of H u ~ d is therefore more precise at larger C values (C ^ 0.3). The 
property of polynomiality is satisfied by construction. Furthermore, we also notice that the 
area spanned by the Bernstein moments, and the value of the missing area calculated from 
Eq. (T4"5l) are in agreement within the theoretical errors. We reiterate that explicit model 
dependence is minimized in the determination of H u ~ d in the ERBL region because of the 
phenomenological constraints obtained from the form factor data, the PDF parametrizations 
for the Non-Singlet (NS) sector, and lattice QCD results. Finally, the results in the ERBL 
region seem to display a visible pattern in the C; t, and X variables. This will allow us 
in principle to construct a specific form for a parametrization driven by phenomenological 
constraints only, that will be present in future work 36]. 

In FigsJT] and [8] we show the results of our parametrization at ( = X for the imaginary 
part of the Bethe-Heitler BH-DVCS interference term at leading order, Eq.([T]), for proton 
and neutron, respectively. The proton results are compared with recent data from Jefferson 
Lab 0|, at XBj — C — 0.36, and Q 2 ~ 2 GeV 2 , while the neutron ones are presented at the 
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kinematics of the forthcoming analysis from Ref. [54| . Notice that Jefferson Lab kinematics 
is expected to be dominated by NS contributions and it is therefore ideal for a comparison 
with the model developed here. We observe an excellent agreement with the data for all 
variants of our parametrization, except for perhaps the largest value in t, that lies, however, 
at the margin of the experimental acceptance. While the proton contribution is governed 
by H, the neutron one is directly sensitive to E, as shown also in our calculation using two 
variants of our fit that differ sensibly in the treatment of the less constrained E distribution. 
We also tested the role of PQCD evolution in this region by showing results at two different 
scales, Q 2 = 2 and 4 GeV 2 , respectively. As expected in the NS sector, PQCD evolution 
turns out to be slow, and it is therefore not to be considered a major concern. 

IV. CONCLUSIONS 

We presented a parametrization of the unpolarized GPDs, H, and E, constructed by 
using a spectator model for the quark-nucleon scattering amplitude in which the forward 
limit is obtained by fitting directly to PDFs in the non singlet sector. This is at variance 
with exisiting phenomenologically constrained models where the forward limit is obtained 
using directly the results from exisiting parametrizations of PDFs, and it allows us to better 
study the role of Regge-type exchanges, that are disengaged, in our case, from the specific 
form used for the parton distributions. The other constraints defining our parametrization 
at zero skewness are provided by the electric and nucleon form factor data. By analyzing 
directly the vertex structure of the scattering amplitude we extended our predictions to the 
non-zero skewness case obtained by modifying the kinematical variables' dependence in our 
expressions for the quark scattering (DGLAP) dominated region. In the X < £ (ERBL) 
region, where a quark- ant iquark pair from the initial proton participates in the scattering 
process, we took into account additional constraints provided by the higher, (^-dependent, 
moments of GPDs. These are in principle available from recent lattice calculations, provided 
reasonable chiral extrapolation procedures are implemented. A reconstruction procedure 
based on Bernstein polynomials, that allows to extract in a largely model independent way 
the GPDs in the elusive ERBL region, is proposed. We showed that with a larger number of 



moments (up to n — 8), a perspective perhaps possible in the near future 55], one could in 
principle obtain a highly accurate description of GPDs in a wide portion of the phase space. 

Differently from the double distribution based approaches that represent a model calcu- 
lation, the proposed method is, for the first time to our knowledge, an attempt to obtain 
a realistic parametrization constrained by phenomenological input, namely by experimental 
data on form factors and PDFs in combination with lattice results. Given the paucity of 
current direct experimental measurements of GPDs, our goal is to provide more stringent 
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and less model dependent predictions that can be useful for both model builders, in order 
to gain insight on the dynamics of GPDs, and for experimenters in planning future DVCS 
and other related hard exclusive scattering experiments. 

Finally, our approach can be compared to recently developed models relying on t- 
dependent moments, namely the dual model, a project to be pursued in future investi- 
gations. In addition, a more detailed treatment of the spectator model used in developing 
our parametrization is on its way, that will allow to extend our approach to the spin de- 
pendent sector where many lattice results exist, and where a rich experimental program is 
being pursued. 
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APPENDIX A: EXTRAPOLATION OF REGGE TERM IN THE DGLAP RE- 
GION AT C + 

At C 7^ 0, the contribution to the parametrization R{^ 2 ) i n Eqs.( l38a] I38bl) . although not 
governed by the kinematical shifts affecting the diquark model kinematics, is modified with 



respect to the ( = case presented in 17| through a ^-dependent kinematical shift in t. The 
shift results from the connection with the formulation of GPDs in impact parameter space 
[53I . where one defines: 

q(x,£,b) = y"d 2 De^ Db if(x,^D) (Al) 

where D, the Fourier conjugate to the impact parameter b, is defined as 

p/ p 

D = — — , (A2 



and: 



4£ 2 M 2 
1 -£ 2 



+ (i - r)D 2 = -t min + (1 - e 2 )D 2 (A3) 



(D reduces to A at ( = [53]). In the overlap representation: 
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By considering gaussian type wave functions, with parameter A, one obtains by Fourier 
transformation, a form for the profile function consistent with Eq. (142 L namely: 



H{x,£,t) oc e" D2/A(a; ' f) (A5) 

where: 



^,0= . a (A6) 



2A 2 1 
1 - x 1 - f 2 

Therefore, for £ ^ one expects the modification with respect to the £ = case of the 
variables in the profile function to be given by: 



^(-t-T^—* (A7) 



1 -x -> (1 -z)(l -£ 2 ) (A8) 
Keeping into account both modifications simultaneously yields the argument of Eq. (l42l) : 



- t -> -t - -^-r = -t + t min (A9) 
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FIG. 1: Left: Amplitude for DVCS at leading order in Q 2 . The light cone components of the 
momenta for the active quarks and nucleons are explicitly written; Right: Time ordered diagrams 
for DVCS: (a) dominant contribution in X > ( region; (b) a qq pair is first produced from the 
nucleon and subsequently interacts with the photons. This process dominates the X < £ region. 
The crossed-terms where two of the particles in the same class are switched, are not shown in the 
figure. 
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FIG. 2: (color online) The dipole masses squared for n = 2, for the isovector magnetic (lower panel) 
and electric (upper panel) contributions obtained by performing fits to the lattice results of [10( 
(full circles). The lines represent the result of our fitting procedure using the ansatz in Eq.([30|), 
The value at the physical pion mass obtained from our fit is also shown (star). 
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FIG. 3: (color on-line) Moments (t) plotted vs. —t, for n = 1,2,3, obtained with the val- 
ues from Table [H compared to the parametrization from Ref . [l7| . For n = 1 (upper panel) the 
parametrization for the form factor from Ref. 17|] is also shown along with the theoretical error 
band. 
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FIG. 4: (color on-line) Comparison between H u ~ d calculated at £ = 0, for different values of t 
(0 < — t < 1 GeV 2 ), using the parametrization in Ref. 13] (full curve), and the reconstructed 
function using the first eight, n = 7, and three, n = 2, Bernstein moments. 
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FIG. 5: (color on-line) Comparison of H u d for different values of £ = 0.18,0.25,0.36,0.53, —t = 
t min = 0.035,0.073,0.18,0.53 GeV 2 (left panel), and -t = -1 GeV 2 (right panel), calculated using 
the parametrization from Eqs. (|38ap . (I39p . (I44p (full curves), and the Bernstein polynomials based 
reconstruction using the first three Bernstein moments. The latter were calculated from the chirally 
extrapolated lattice QCD moments (Figj3]), at the scale Q 2 = 4 GeV 2 . The horizontal error bars 
are obtained by evaluating the dispersion, Eq. (|34p . The error bars on the moments are mostly not 
visible on the plot. The shaded area represents the "missing area" in the ERBL region, obtained 
imposing the polynomiality condition (see text). 
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FIG. 6: (color on-line) H u ~ d in the ERBL region obtained with the method described in the text, 
for different values of ( = 0.18, 0.25, 0.36, 0.53, and 0.10 < -t < 1 GeV 2 . 
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FIG. 7: (color on-line) The imaginary part of the Bethe-Heitler BH-DVCS interference term at 
leading order, C J (F), Eq.([T]) for a proton target. Experimental data are from Ref. 26] at XBj = 0.36 
and Q 2 = 1.5, 1.9, 2.3 GeV 2 . The full curve shows our prediction, including the theoretical error, 
evolved to Q 2 = 2 GeV 2 . The dot-dashed curve, obtained at Q 2 = 4 GeV 2 , shows that the effect of 
evolution is relatively small. The other curves represent variations of the parametrization obtained 
respectively by disregarding the contribution from E, in Eq.([T]), by disregardin g th e kinematical 
shift in the Regge term described in the text, and by using E from Set II of Ref. 171 ]. 



31 



f 0.4 
0.2 





-0.2 
-0.4 
-0.6 
-0.8 




— AHLT, Set I 
AHLT, Set II 



-0.5 -0.45 -0.4 -0.35 -0.3 -0.25 -0.2 -0.15 

t (GeV 2 ) 

FIG. 8: (color on-line) The same as FigJT] for a neutron target. Calculations were performed in 
the kinematical range of Ref. 54], namely at Q 2 ~ 2 GeV 2 and xbj = 0.36. The two curves are 



variations of the parametrization using E from Set I and Set II of Ref, [171] . respectively 
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